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Abstract: The simplest version of a class of toy models for QCD is presented. It 



^ , is a Lipkin-type model, for the quark-antiquark sector, and, for the gluon sector, 

^ ' gluon pairs with spin zero are treated as elementary bosons. The model restricts 

to mesons with spin zero and to few baryonic states. The corresponding energy 
^ ' spectrum is discussed. We show that ground state correlations are essential to 

■ describe physical properties of the spectrum at low energies. Phase transitions 

are described in an effective manner, by using coherent states. The appearance 
of a Goldstone boson for large values of the interaction strength is discussed, as 
related to a collective state. The formalism is extended to consider finite tem- 
peratures. The partition function is calculated, in an approximate way, showing 
the convenience of the use of coherent states. The energy density, heat capacity 
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and transitions from the hadronic phase to the quark-gluon plasma are calculated. 
PACS: 12.90+b, 21.90.+f 

1 Introduction 

QCD is considered to be the theory of the strong interactions. It is well 
understood at high energies. At low energies, the QCD coupling constant 
becomes too large to apply perturbation theory. Lattice gauge calculations 
Ijl]] may describe the non-perturbative QCD regime, instead. Yet, problems 
like finite size effects and fermion doubling still persist. Although some 
advances have been made, only the lowest states for a given spin and par- 
ity can be calculated. The observed sequence of levels cannot be explained 
by lattice gauge calculations and alternative methods have to be develop 
to explain the ordering |Q. Many effective models have achieved some suc- 
cess in describing the low energy regime of QCD 0, ^. These models have 
in common that only quarks and antiquarks are taken into account in the 
fermionic sector, while effective gluon potentials or states with a fix number 
of gluons are considered. In the real world hadrons are built by quarks, an- 
tiquarks and gluons |^. The interactions between these degrees of freedom, 
in consequence, may play an essential role in order to understand QCD at 
low energies. 

Concerning QCD effects at finite temperature, i.e. the investigation of 
the quark-gluon plasma (QGP) 0, it exists an intense effort, mainly focussed 
on computational aspects of the problem 

In this paper, we propose a toy model for QCD which: i) is amenable 
for an analytical treatment, except for a numerical matrix diagonalization, 
ii) may describe the meson spectrum for flavor (0, 0)-spin and partly the 
baryon spectrum, iii) it is able to describe phase transitions at T = 0, as 
a function of coupling parameters, iv) can describe some characteristics of 
the transition from the hadron gas to the QGP , and v) can be used to test 
microscopic many-body techniques, intended to describe realistic scenarios 
of QCD which cannot be accessed by other methods. 

The model is meant to mock up the basis features of non-perturbative 
QCD, in a similar way as some schematic models do in the nuclear many- 
body problem |10]. We should stress that our toy model does not result 
from a field theory, rather it represents the interactions between effective 
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degrees of freedom, instead. The basic ideas, procedures and methods can 
be discussed already for the simplest versions of the toy model in a very 
transparent way. At a certain point in our discussion we shall concentrate 
on the structure of an extended version of the model, which can be treated 
analytically. 

The paper is organized as follows. In section 2 the model is introduced, 
the energy spectrum is calculated and its structure is discussed in terms 
of the elementary degrees of freedom. Coherent states are introduced to 
determine the occurrence of phase transition, induced by variations of the 
strength of the interactions. In section 3 we discuss finite temperature effects, 
by introducing temperature and by calculating the grand canonical partition 
function. Conclusions are drawn in section 4. 



2 The Toy Model: Zero Temperature Case. 

The fermion sector of the theory is described by the action of the operators 
which create (annihilate) quarks, with effective masses Uf. Schematically, it 
corresponds to the situation represented in Figure 1, where two levels, with 
energy ioj/ and degeneracy 2f2 are represented as valence space [0]. The 



degeneracy of each level, is given by the product of the number of colors (rzc), 
spin ijis), flavors ijifi) and all other possible degrees of freedom {ricoi), like 
orbital quantum numbers, etc.. For temperature T = and no interaction 
the lower level is filled by fermions. The creation (annihilation) operators of 
these fermions are c^^^^y q^^^. (c'^(^''^)/°'*)^ in co- and contra- variant notation for 
the indices. The symbol (1,0)/ refers to the flavor part, where (1,0) is the 
SU(3)-flavor notation and / is a short hand notation for the hypercharge Y , 
the isospin T and its third component T^. The index a represents the two spin 
components ±i, the index i = 1 or 2, stands for the upper or lower level and 
the index a represents all remaining degrees of freedom, which are at least 3 
if the color degree of freedom is considered. Lowering and raising the indices 
of the operators introduces a phase, which depends on the convention used 



13| , and a change of the indices to their conjugate values, i.e., the quantum 
numbers (1, Q)YT%a change to (0, 1) -YT -T^- o. 

The quark and antiquark creation and annihilation operators are given 
in terms of the operators c and 
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which corresponds to the Dirac picture of particles and antiparticles: quarks 
are described by fermions in the upper level and antiquarks by holes in the 
lower level. 

The gluon sector of the model space is described by bosons which rep- 
resent pair of gluons coupled to spin zero. The energy of a boson state is 
fixed at the value Ub and the state is created (annihilated) by the action of 
a boson creation (annihilation) operator b on the vacuum. 

The quark-antiquark pairs of the model are given by 



a a 

»-^/i<7i2 - -0/1(71 - 2^ c^/i<7i2<= -2^ "a/101 a 

a a 

W20-21 _ „t „a/20-2l _ t „a/20-2 

*-^/ioil - Z^'^a/ioilC -Z^"a/i(7i" 

a a 

/-«/20-22 _ t „a/2022 _ J jta/202 /ON 

*-'/l(Ji2 — 2^^aha^2^ — 2^ "'afi'^i"- ■ K"^) 

a a 

The first two equations describe the creation and annihilation of quark- 
antiquark pairs. The pairs can be coupled to definite flavor (A, A) = (0, 0) or 
(1, 1) and spin S" = or 1. We shall write, in this coupling scheme, -B|;j^ 
where / is the flavor, S is the spin and M is the spin-projection. The 
operators B(^x,\)f,SM annihilate the vacuum |0 >, which is the configuration 
where the lower state is completely filled and the upper one is empty. The 
operators in (H) form a U(12) algebra. To simplify the discussion, we shall 
restrict to a sub-algebra given by the pair operators coupled to flavor singlet 
((0,0) in the SU(3) notation lH), 



5+ — a/6-B|oo)ooo 
S- = a/6B(o,o)o,oo 

5o = Uf-Q , (3) 

with rif = l!l9±!]sl^ where Uq is the number operator for quarks and riq the 
number operator for the antiquarks. They form a SU(2) algebra. The model 
is equivalent to the Lipkin model P], familiar in nuclear physics, with the 
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difference that the operators are given here by the combination of quark- 
antiquark pairs . The addition of an extra boson level was also discussed in 
Ref. |]10[ and it is related to pion effects in nuclei. 



As Hamiltonians we shall consider two different types: 
Hi = 2ufSo + uJbnb + Vi{Slb + b^S'i) 

Hji = 2ufSo + ujtnf, + Vi:{S+ + S-)'':{b^ + b) , (4) 

with So = Uf — Q and the double dots indicate normal ordering. 

The Hamiltonian Hj exhibits a useful symmetry, i.e. it commutes with 
the operator 



^ + . (5) 



This version of the Hamiltonian is similar to the one given in Ref. |TD[, 
except that in the SU(2) operators appear linearly while in (^) they 
appear quadratically. Because a quark-antiquark pair has negative parity, 
the operators and should appear quadratically to conserve parity. 

The states of the model space, belonging to Hj, are SU(2) states with 
the additional ordering given by the eigenvalues of P. The vacuum state is 
defined via 6|0 >= S-\0 >= 0. Because the number of fermion pairs Uf is 
limited by 2Q the range of ni, for a fixed value of P is also limited. Therefore, 
the matrix representations of Hj are finite. A large eigenvalue of P implies 
that the corresponding configuration has many gluons. In Figure 2 we show 
the energy of the lowest state for a given strength of the interaction, as a 
function of the eigenvalue of P. For zero interaction, the energy increases 
monotonically. For large values of Vi it appears a minimum with a large value 
of P. This implies that the physical ground state should be a correlated one, 
that is to say that the physical ground state will, likely, be a state with a 
large number of gluon pairs. Concerning the dependence upon Vi, which 
is the strength of the interaction which couples pairs of gluons with pairs 
of quark-antiquark pairs, the curves of Figure 2 show minima, which are 
different from the perturbative vacuum, for values of Vi > 0.035 GeV. For 
larger values of the interaction strength Vi, the lowest state is the one with a 
large eigenvalue P, indicating a gluon dominated vacuum (for the transitional 
region, quarks- ant iquark and gluon pairs could appear in the vacuum with 
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comparable weights). In Figure 3, the energy spectrum of Hj for positive 
parity states is displayed. The parameters used in the calculations are: i7 = 9 
(i.e., nc = 3 and Ufi = 3), = | GeV and Ub = 1.6 GeV. For Vi = GeV 
we obtain the first state at | GeV, corresponding to two quark- ant iquark 
pairs. The next state is the glueball at 1.6 GeV. When the interaction is 
turned on, the energy changes until it reaches a "critical" or "transitional" 
point at V{ = 0.035 GeV. There, a level crossing occurs and the lowest state, 
for higher values of Vi, has both quark- antiquark and gluon pairs. Beyond 
the transitional point, the density of levels increases. This effect is known 
from nuclear physics, where the transition from a spherical nucleus to a 
deformed one is accompanied by a significant increase of the density of levels 



at very low energies |TH]. For values of the interaction larger than the 
ground-state expectation-value of the number of quark-antiquark and gluon 
pairs increases. For large interaction strength the number nj approaches 
a constant value, reflecting the Pauli principle, i.e. only a certain number 
of quarks can occupy the higher level. This behavior will be discussed in 
detail for the case Hjj. The relatively high density of states at low energy, 
shown in Figure 3, may not be very realistic. However, the model predicts 
the appearance of some states at very low energy. This would correspond 
to a pion-like structure, which is also indicative of a collective nature. This 
problem does not appear in the case of the Hamiltonian Hjj, which does not 
commute with P. The Hamiltonian Hjj has to be diagonalized in the whole 
space, which is infinite dimensional. The diagonalization can be performed 
numerically by introducing a variable cutoff in the number of bosons. We 
have adopted, as a criterium for convergence, the stability of the low-energy 
sector of the spectrum as a function of the cut-off. The Hamiltonian Hji 
contains all terms which are required by symmetry, i.e. it includes a term of 
the form S^b, describing the annihilation of a gluon pair and the creation 
of two quark-antiquark pairs, and also a term S'^b\ describing ground state 
correlations. The scattering term S^S-{b^ + b) appears with a factor of 
2 because fermion lines can be exchanged. Because of the symmetry in 
permuting the lines all interactions should have the same (or at least similar) 
coupling constant, justifying the use of only one interaction parameter, Vi. 
Figure 4 shows the dependence of the energy, of positive (Case (a)) and 
negative (Case (b)) parity states, on the strength Vi, for the values of tuj and 
Uh given above. The values are referred to the positive-parity ground-state. 
Contrary to the case of the Hamiltonian Hj, the Hamiltonian Hn does not 
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show a dense spectrum beyond the phase transition point. In the transitional 
region several avoided crossings occur and the spectrum is richer there, than 
outside that region. Concerning the behavior of the low-energy part of the 
spectrum, it shows, after the transition point, a negative parity state which 
is degenerate with the positive parity ground state. This state, a Goldstone 
boson, can be interpreted as a collective, pion-like, state whose structure can 



be understood in the framework of coherent states |]18| , as it will be discussed 
later on. 

Figure 5 shows the difference of the expectation value of the number of 
gluon- and fermion-pairs, as a function of the strength, Vi, of the interaction, 
for the case of the Hamiltonian Hjj. As it is seen from the figure, the results 
can be interpreted in terms of equal population of fermion and gluon pairs 
{Vi < 0.008 GeV), fermionic dominance (0.008 < Vi < 0.015 GeV), and 
gluonic dominance {Vi > 0.015 GeV). 

The features of the spectrum, and of the ground state occupation num- 
bers, can be understood in terms of a "phase transition" at zero temperature, 
that is to say, in terms of a change in the correlations induced by the Hamil- 
tonian. A convenient way, to represent this effect, is to introduce a set of 
states (coherent states) and to define an order parameter (the ground state 
occupation number associated to a given degree of freedom). The calcula- 
tion of the expectation value of the Hamiltonian, in the basis of coherent 
states, and its variation with respect to the order parameter yields the possi- 
ble "phases" of the system, as extremes of the minimization procedure. The 
set of coherent states, which we have adopted in our calculations, is defined 
by PI 



\z > = \Zf > \zi, > 



,2 , t . 



\zb > = e-^e^"^ |0 >b , (6) 

where |0 >/ and |0 >h are the fermion and the boson vacuum and |0 > is the 
product vacuum state (|0 >/ ®\0 >b). The power Q, of the fermion normal- 
ization factor, indicates that we are working in the SU(2)-spin-representation 
which includes, as the lowest state in energy, the fully occupied level at —Uf 
(see Figure 1). 
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By defining the complex order parameters zj = pfC^'^f and Zb = Pb^^'^'', 
the expectation value of the Hamiltonian Hjj (|^) is given by 

<Hii> = -2Quf + Ubpl 

+ V,i^ + + j2pbCos{ct>bm) 

and it can be regarded as a classical potential in the parametric space of 
the amplitudes p and phases 0. It shows for small coupling constants Vi 
a minimum at \zf\ = pj = 0, which represents small departures from a 
dominant harmonic (quadratic) potential, and a deformed minimum for a 
sufficiently large value of Vi and a given combination of 0/ and (pb- The 
factor 2(pf, makes the potential invariant under the change 0/ — ^ 0/ + it. 
Thus, when the difference in energy of the two minima with respect to the 
barrier between them is sufficiently large, there exist two degenerate states, 
one with positive and the other with negative parity, which minimize the 
expectation value of the Hamiltonian. 

The pseudo-scalar particle at zero energy can not be identified with the 
pion, because the pion belongs to a flavor (1, 1) (octet) representation of 
the flavor group SUj(3). However, the fact that the model has a Goldstone 
boson gives some hope that a generalized model with open spin and flavor, as 
indicated at the beginning, may also exhibit a low lying negative parity state 
which can be identifled with the pion. This conjecture is justifled because 
an equivalent interaction, as in Eq. (^, for flavor (1, 1) spin will exhibit a 
similar behavior with respect to a coherent state which includes pairs with 
flavor (1, 1) spin 0. 

In our model the basis for the baryonic states is given by 

\q'iqqnX,p)f,SM> ~ (5+)"|g^ (A/.)/, > . (8) 

the index (A,/i)/ refers to the flavor, which is (1, 1) for the octet, etc. The 
indicates that the state to the right, on which the operator acts, is a 
pure three-quark state. The three-quarks state satisfles 

S.\qMXp)f,SM> = . (9) 
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That annihilates the three-quark state holds because the quark-antiquark 
pair operator contains an antiquark annihilation operator which anti-commutes 
with the quark creation operators of the state on the right and annihilates 
finally the correlated vacuum. The basis for the mesons is obtained using 
(A,/i) = (0,0) and S' = M = 0. Note that the Hamihonian Hu (and the 
same is true for Hj) does not distinguish between different flavors, and, 
therefore, the flavor- (1,1) and flavor-(3,0) baryons are degenerate. Part 
of the degeneracy can be removed, by introducing terms depending on the 
hypercharge and the isospin. In order to remove completely the degeneracy 
between the (1, 1) and the (3,0) flavor conflgurations one has to include, of 
course, flavor-depending interactions, as it will be done in a more general 
formulation of the present toy model |TI[ . 

The problem for the baryons is completely analogous. Due to the fact 
that three quarks minimally occupy the higher level of the fermion model 
space, the effective degeneracy, i.e. the number of conflgurations available 
to excite quarks from the lower level, is 21] — 3 (since the total number of 
available states is 2Q). The factor 3 is a consequence of the Pauli-blocking. 

In Figure 6, the spectrum of the lowest baryonic states, as a function of 
the coupling strength Vi and referred to the lowest positive-parity mesonic 
state, is shown. After the transition point, the energy of the baryonic states 
increases. To obtain states below 1 GeV, one has to reduce the effective 
quark mass. Since we are interested in the trends exhibited by the spec- 
trum, we shall not flt it to physical masses. The transition point is slightly 
shifted to higher values of the coupling constant, due to the lower value of 
the degeneracy, otherwise the behavior observed for the baryon spectrum is 
similar to the meson case. The relatively delayed onset of a transition in the 
baryon spectra produces a region where the physical vacuum may contain 
a pair condensate while the baryonic states may still be built as pure three 
quark objects. 

In Figure 7 the content of quark-antiquark and gluon pairs in the flrst- 
excited barionic state is shown. The results show a large contribution from 
gluon pairs to the baryonic states after the transition point. This indicate 
that in a more realistic model and for sufficiently large interaction strength, 
some baryon states, like the proton, may contain a sizeable contribution due 
to sea quarks and gluon pairs. 

After the analysis of the properties of the low lying spectrum of the two 
Hamiltonians (^, we can proceed to discuss finite temperature effects, which 
may be relevant for the description of the transition from the hadronic phase 
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to the QGP. 



3 The toy model: Finite Temperature Case. 

In this section we shall present the results corresponding to the finite tem- 
perature case. We shall show the main steps related to the calculation of the 
partition function, which has been performed by extending the techniques 



discussed in \16 



In the limit Vl = 0, the fermionic sector of the Hamiltonians Hj and Hji 
reduces to the free Lipkin model, which consists of two levels, with ener- 
gies ±ujf and a degeneracy (2fi). Allowed configurations are specified by all 
possible arrays of particles and holes, in both levels, and their degeneracies. 
Thus, a certain configuration can be specified by listing the number of occu- 
pied (empty) states in the upper (lower) level. The operators which create 
(annihilate) these fermion pairs obey a pseudo-spin block-algebra, for spin 
i. There are in total 2Q building blocks. If i>i denotes the number of blocks 
where both levels are occupied, 1^2 is the number of blocks where both levels 
are empty and 2r is the number of blocks where either one of the levels is 



occupied, the partition function can be written as |17] 



m = i: p^izft'/.-. . (10) 

where 

with J = T — k, P = ^, and T is the temperature in units of GeV. The 
states \z >= \z >f >i, are the normalized coherent states \TE\. The 
factor gl is the multiplicity of the configuration with pseudo-spin J = r — k. 
For each value of J one should define a coherent state \zfj >. The coherent 
state, used in the previous section (the T = case), corresponds to the value 
J = Q. The chemical potential /i multiplies the operator N, which gives 
the total number of particles in the lower and upper level. The partition 
function (|l^) does not conserve fiavor, color and spin. A multiple projection, 
to restore these symmetries, can be carried out, in principle, although it is 
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a very involved procedure [0. We assume that the volume occupied by the 
QGP is large enough so that a projection is not needed, though we shall 
restrict to a sub- volume for the thermodynamic description. 

If interactions are added to the Hamiltonian, the partition function cannot 
be obtained analytically, in general. If the value of Q is not too large, the 
partition function (p!OD can be obtained numerically. 

We have diagonalized the Hamiltonian Hj and obtained a set of eigen- 
values for each value of J. The parameters of the Hamiltonian were fixed 
at "l/i = GeV and 0.04 GeV, = 9, cj/ = | GeV, and Ub = 1.6 GeV. 
The partition function ( [I0| ) was calculated for temperature T < 0.5 GeV. 
Since the value of Q is not very large, one may ask if, at high temperature, 
more configurations (i.e. larger values of fl) should be included. We have 
checked, numerically, this effect upon the partition function. Figure 8 shows 
the value of (pUD, for each value of J < fl. As seen from this figure the 
contributions to the partition function reach a maximum for a certain value 
of J, which is smaller than Q. This result justifies the approximation with 
J >> 1. First and second derivatives of the partition function (the internal 
energy and the heat capacity) are shown in Figures 9 and 10, respectively. 
The heat capacity, for Vi = GeV, shows the Schottky bump typical 
of a two level system. The shape of the curve remains unchanged when the 
interaction is switched on in Hj. The increase of the interaction strength 
Vi produces a sharper peak in the curve, indicating a possible second order 
phase-transition. The results shown in Figures 9 and 10 have been obtained 
with the Hamiltonian Hj. Although we are not showing it in the Figures 
9 and 10, the same thermal behavior of the heat capacity is obtained with 
the Hamiltonian Hu. However, the values of (T and Vi), where the onset 
of the phase-transition is produced, are different for the partition functions 
corresponding to Hj and Hjj. 

The above discussed results, which have been obtained by performing a 
numerical diagonalization are indeed the exact results of the model. Poten- 
tially, they exhibit the desirable thermodynamical features of QCD. We shall 
take these results as reference values for an approximate calculation. The 
obvious motivation for such approximate treatment is the generalization to 
larger values of the model parameters. 

We shall discuss first the treatment of the fermionic sector, i.e. the inte- 
gration on the coherent states | ^^(/j) >• For large values of J, the overlap 
< Za\zf3 > is a decreasing function of \pa — Ppl- Also in this limit, it is a 
strongly oscillating function of the phase difference <pa — <t>i3- Due to this 
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oscillation we can write, both for and 



< Za\e '^^\zp > < Za\e '^^{Za >< Zo\zp > , (12) 

where the overlap in the above equation takes into account the off diagonal 
behavior for a ^ (3. The expectation value < Za\e~^^ \za > can be expressed 
as an integral product of the form 

TV J (1+|2^|2)2) ^ ^a\C \Za ^ — 

y-oo (-/3)" (2J+1)" r d^z^^...d^z^^_^ ^ ^ I I I I 

Z^n=0 n! vr" J n)!~"^{l+|2-yj. p)2 ^7ol-*-*l^7i I I 72 ^ 

••• < 2:^„-il^fk7o > ) (13) 

by inserting n-times the unit operator . The Hamiltonian H is expressed 
into two terms Hq + H', where Hq is the non-interacting part and H' 
contains the interactions. The k-th term of the expansion has k matrix 
elements of H' and (n — k) factors with Hq. By performing a rearrangement 
the contribution of the k-th partition to the integral is written 

g {-pr (2 J + 1)' f (Pz^cPzp ( n 



n=k 



J (1+|^,|2)2(1 + |;,^|2)2 ^ k 



<z^\{Hr\zp><zp\{Hor-'\z^> . (14) 



This result illustrates the binomial character of the k-th partition. By a 
re-exponentiation one can write 

TT^ y (1 + |^„|2)2(1 + |^^|2)2 < X 1^" > \^^) 

The integral / dz^^ J dz^, can be replaced by a single integral, by using equa- 
tion ([T^ ) and by performing the integration on the variable z^. The validity 
of this approximation is restricted to large values of J, whose validity was 
shown above. After these approximations (|13D reads 

r d^z^d^Z0 I ^ lp-/3/fV > 

7r2 J (l+|2„|2)2(i+|2^|2)2 ^ Za\e \Zf3 >^ Zf3\e \Za > 
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Since Hq = 2LjfSo, the matrix element < Za\e ^-^"\za > is readily 
calculated and the result is 



(1 + Izar)'^^ 



As a check on the consistency of these approximations, it is verified that 
by setting H' = and from the above equation, the integration of (|lT]) gives 
a result which is identical to the one of . 



Concerning the matrix element of the interaction H', we have adopted 
the following approximation 

< Zo^le-f"^' \za > ~ e-^<^"l-^'l"'^> , (18) 

which is valid when the temperature is high and/or the interaction coupling 
constant is small. It corresponds to the factorization < Za\{H')"'\za > ~ 
(< Za\H'\za >)", a result which can be reproduced, by assuming (|l^) and 
inserting n-th unit operators between factors H' . As said before, the applica- 
bility of the procedure is limited to relatively small values of the interaction 
strength Vi and relatively large values of the temperature. 

Let us now turn the attention to the bosonic degrees of freedom. The 
exponent of Eq. ([18|) is a linear combination of the boson operators and 
b. The coefficients are given by the expectation values of or S+S- (see 
appendix A), which are functions of the complex variables Za and z^. The 
normalized boson coherent state is given by |]TB| 



\zb > = e-^e'^b |0 > . (19) 



The volume element of the complex integral is ^ (in Ref. [M the coherent 
states are not normalized and therefore the volume element has an extra 
factor e"'^*' ). The calculation involves an integral of the type 

- / d^Zh < Z6|e-^"'^''+"i^'+'^2^|z6 > . (20) 

IT J 

The coefficients a^, which depend on the expectation value of powers of the 
SU(2) generators and which are proportional to the interaction strength Vi, 
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have a common value for the model Hamiltonian but they differ for the 
case of Hi. 

In order to evaluate the expectation value which appears in (pO]), the 
Baker-Campbell-Hausdorff formula is applied 0, 0. Notice that one can- 
not apply the approximations described for the fermion sector, because the 
overlaps < zi,\z[ > show a broader dependence on — p[ and — (f)[. The 
result is 



with 



K 

and where 



66 



€f + 6 



2l3Va 



(l-e-'3"(.) 



, 6(//) 



2 + /3a;fe 
(1 - e-P'^b) 



(22) 



K = Vi< z^\Sl.\z^> 

Vc = Vi< z^\ : (5+ + ■.\z^> . (23) 

The indices I and // refer to the two different model Hamiltonians. 

Because the exponential function which contains the operator b is acting 
on the coherent state \zb >, the operation is well defined and it gives a 
factor e^^^^ The same holds for the exponential function which contains the 
operator and it gives e^^^b. The expectation value of the operator which 
contains reads 

< Zb\e^''^'\zb > = e-^^'^'e^^'^'^'' . (24) 
Finally, the integration over the complex variable Zb yields 
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- [d'z< z,|e^^^'e«^^^e^^^|2, > = . (25) 

7T J 1 — 6^2 

In the following we shall construct the final expression of the partition func- 
tion, after the above introduced approximations. We shall restrict to the 
case of Hj, in order to compare the results of the approximations with the 
numerical (exact) results. A similar analysis can be performed for the case 
of Hjj The remaining integration on the complex variable Za, of the 

fermion part, leads to 



(2J+1) r dx , ^''7' , (1 + .e-^^^fre'''''''^^^'^ ,(26) 



with 



F{(i) = 2 — ^ ^ ' 



uj (2 + (3ub){e(^''^ - 1)(1 - e-(^^b) 

after substituting x = p^. F{[3) is a smooth function of (3, which approaches 
the limit ^ for large (5. The integral, (^), is a function of the pseudo-spin 
J and its argument is the product of exponential functions of positive and 
negative functions of J. That the integral can be represented by the integrand 
at J = Jo, where Jq is the value that maximizes it, can be easily seen by 
setting Vi = 0. 



Let us call g{J) the ratio between the integral (26) with Vi ^ and 
the same integral with Vi = 0, and Jq the value of J which maximizes the 
integral. We can write the Taylor expansion of g{J) as 




9iJ) = (^-^o)"- (27) 

By using the identity 



rx^^ = (28) 
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with X = e'^'^f, the expansion in terms of {J — Jq) of g{J) transforms into 
an expansion in the derivatives Thus, the sum (pT]) is of the form 

J2n D^{Joi ^)Zo, where Zq is the partition function for the non-interacting 
case fl^ and D is the differential operator (y'Sc ~ "^o)- choose Jq such 
that the ~ Jo)Zo = 0. The value of Jq, obtained in this way, is given 

by 



Jo ^ n 



X^Y \ ,^ 1 



(x + y)(i + xy) / ' 



where Y is the fugacity [|T^. This expression approaches Jq = f2 for /3 — > oo 
and it gives small values of Jq for /3 ^ 0. After these approximations the 
partition function can be written as 



Z — ZfZbZint , (29) 

where Zf and Zi, are the partition functions for free fermions and free bosons, 
respectively, and Zint is the contribution due to the interactions. For (3^0 
Zint ^ 1 by construction, and the high temperature limit is automatically 
obtained. 

In what follows we choose the same values of interaction strength as done 
above when the partition function was calculated numerically. 

In Figures 9 and 10 we show the internal energy and the heat capacity re- 
spectively, as a function of the temperature, obtained with the approximate 
partition function (^). As it can be seen from the curves, the approximation 
works reasonable well for high temperatures, T > 0.15 GeV. At low tempera- 
tures the curves, representing the approximated values, deviate significantly 
from the exact results. 

The behavior of the curves of Figure 9, illustrates the agreement between 
the approximated and the exacts results. It is valid only for the high tem- 
perature region. At low T, the main contribution to the partition function 
comes from the ground state. The lowest energy is the one of the configu- 
ration where all states in the lower level are occupied while the ones in the 
upper level are empty. This corresponds to the case (2r) = 2Q, ui = U2 = 0. 
The partition function reduces to one term, given by the integral ( pB]) with 

J = fi, and the energy goes to — as T goes to zero. This value, is of 
the order of —13 GeV if one uses the already given parameters (Vi = 0.04 
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GeV, f2 = 9 and Uf, = 1.6 GeV). The exact calculation yields —2 GeV. This 
deviation is caused by the approximation ([T8|). 

Note that, in a more realistic context, the transition to the QGP is be- 
lieved to take place around T = 0.165 GeV or T = 0.270 GeV For these 
temperatures our approximate results are acceptable. Also, for these tem- 
peratures there is still a sizeable difference between the results obtained with 
interactions and without interactions. It shows also that the interactions 
cannot be neglected in the high temperature regime. 

Finally, we like to discuss the dependence of T on the chemical potential 
fi, when the pressure of the system is equal to the bag pressure B = 0.145GeV 



21| . To find such a dependence one has to introduce a volume. We consider 
an elementary volume given by the size of a hadron Ifm^), as done in 
0. The pressure is then given by the ratio of the internal energy and this 
volume. The result is depicted in Figure 11. Without interactions (Vi = 0) we 
reproduce the results of Ref. ||21|. By turning on the interaction (Vi = 0.04 
GeV), the chemical potential increases. This effect shows that the correlated 
vacuum state is dominated by gluon pairs and it also has contributions from 
quark- ant iquark pairs. 



4 Conclusions 

In this work we have introduced the essentials of a toy model for QCD. The 
model consists of two levels with energies ±u;/, which describe the fermion 
degrees of freedom, and gluons are introduced via a level of positive energy 
which can be filled by gluon pairs with spin zero. The gluon pairs are treated 
as bosons. Two different Hamiltonians are discussed. The first one, (Hj), 
commutes with the symmetry operator P, a property which allows us to cal- 
culate the energy spectrum easily. However, this Hamiltonian contains only 
a certain type of interaction terms which do not include ground state corre- 
lations. The second Hamiltonian (Hji) does contain terms which produce 
ground state correlations. Due to the symmetry of the vertices, respect to 
the exchange of fermion and boson lines, all terms entering in Hjj have the 
same interaction strength. We have shown that the corresponding spectrum 
exhibits a phase transition, depending on the interaction. For small values 
of the interaction Vi, fermion pairs and gluon pairs equally populate the 
ground state, for intermediate values of Vi the physical vacuum is described 
by quark-antiquark pairs while for larger values of Vi gluon pairs dominate. 



17 



In the gluon dominated phase the spectrum has a degeneracy of the ground 
state, given by one positive- and one negative-parity state. This property 
gives some hope that a more general version of the model, with open flavor 
and spin channels, may show a Goldstone boson ( in the flavor octet (1,1) 
SU(3) notation), if the strength of the corresponding interaction is large 
enough. The appearance of the gluon and quark- ant iquark condensate and 
of the Goldstone boson, may be easily described by using coherent states, as 
we have shown for the present version of the model. The same approach may 
be useful for the case of more general models of non-perturbative QCD. 

Baryons were also considered in the model (see Eq.(8)) and due to Pauli- 
blocking effect the effective degeneracy Q decreases. Also, the transition 
point to a condensate of pairs of quark- ant iquarks and gluons is shifted to 
larger values of the interaction strength. There is a regime where mesons are 
already in a condensate phase while the baryons can still be treated as three 
quark systems. We have investigated finite temperature effects, by construct- 
ing the partition function of the model, both exactly and approximately. We 
have shown that the use of coherent states makes it possible to introduce ap- 
proximations in a controlled way. The results, for the internal energy, heat 
capacity and the equation of state of the system at the bag pressure, are in 
agreement with previous calculations. To summarize, the model is able to 
describe characteristic features of QCD at low and at high temperature. This 
gives some hope that, in a more general version of it with open flavor and 
color, it may describe the hadron spectrum at low energy and the transition 
to the QGP, as well. 



Appendix A 

We want to calculate 
and S+S- using coherent states. We adopted, as a suitable representation. 



We want to calculate the expectation values of the operators So, (>S'i + S'i 



the normalized coherent states \18 



I Zf >= ^^6^-^^+ \ J,-J > 

where | J, — J > is the eigenstate of the pseudo-spin algebra with the absolute 
value of the spin J and its projection M = —J. 
For the operator Sq we have 
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Now, since [S-, So] = S- and 



+ i[A,0] + i[A, [A,0]] + 



(30) 



we arrive at 



e'f^-Soe'f^+ = {So + z}S^)e'*f^- e'f^^ 



This yields 



<z,\So\z, >= <J,-J\ (So + z}S_)e''f^-e''^* | J, -J >= 



__i_(_ J + ,;^) ,-;S-e'/«. I J, -/ > 

* s s 

To calculate the normalization, < J, — J | e^^ -e^/'-'+ | J, — J >, we write 
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2J 



(l + |z^|)-^|z^>=e^/'^+| J,-J>=^ 



n=0 



(2J)! 
\\ {2J~n)\ 



J,n > 



leading to <J,-J\ e'f^~e'f^+ \ J, -J >= (1 + \zf\y-^. 
With this result, the expectation value of So reads 



< Zf \ So \ Zf > 
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Next we are going to calculate the expectation value of the operator 5^ 



< ^/ I Si I z, >= ^Y^^^ < J, -J I (e^f^-Sle^fS-) I J, - J 



1 



(1 + \zf\^r^ dz} 



<J,-J\ {e'fS-e'i'S+) I J,-J>^ 



Therefore 



< Zf I + I Zf >- jYTJ^^ 

To calculate the expectation value of S+S- we use the identity S-^-S- 
S'^ — Sl + So, and from the above obtained results we find 

<z,\Sl\z,>^^^^^i-J^z^,^)^^^ 



, _ AJ'\zf\' 2J\zf\' 2J{2J-l)\zf\' ^ 2J{2J-l)\zf\' 
and 



< Zf I S+S- I Zf >=< Zf I (5^2 - S'^ + 5o) I ^/ > 



2J(2J-l)|z,P J{l-\zf\') _ 2J\zf\H2J+\zf\') 

Finally, in the limit of large J expressions like (2 J — k) are approximated 
by 2 J. 
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FIGURE CAPTIONS: 



> 

CD 



-CO, 



Figure 1: Schematic representation of the model space. The fermion levels 
are indicated by their energies ±<^/. The gluon-pairs are represented by the 
level at the energy u>b. 
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V, =0.000 GeV 
Vi = 0.020 GeV 
V, = 0.035 GeV 
V,= 0.050 GeV 



Figure 2: The energy of the ground state, £"0, in units of GeV, as a function 
of P, for several values of the parameter Vi. The values correspond to the 
calculations performed with the Hamiltonian Hj. Note the occurrence of a 
nontrivial minima when Vi > 0.035 GeV. 
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Figure 3: The spectrum of the model Hamiltonian Hj, for positive parity 
states, as a function of the couphng parameter Vi. Note the crossing of 
excited states with the perturbative ground state at about Vi ~ 0.035 GeV. 
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Figure 4: The energy spectrum of Hamiltonian Hn for positive (Case (a)) 
and negative parity (Case (b)) states, as a function of the couphng strength 
Vi- The on set of the phase transition takes place for lower values of Vi, as 
compared with the results corresponding to the Hamiltonian Hj (see Figure 
3). 
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Figure 5: The difference between tlie vacuum expectation values of the 
number of gluon pairs rif, and fermion pairs Uf is shown, as a function of Vi, 
and for the case of the Hamiltonian J///. 
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Figure 6: The spectrum of the lowest baryonic states as a function of the 
interaction strength. The values of the excitation energies are taken with 
respect to the ground state of the mesonic sector. The Hamiltonian Hfj was 
used. 



28 



0.005 0.01 0.015 0.02 

Vi [GeV] 

Figure 7: The content of quark- ant iquark (nj) and gluon (n^) pairs in the 
first excited state of the barionic spectrum, as a function of the interaction 
strength Vi. The Hamiltonian Hjj was used. 
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Figure 8: Contribution to the partition function, Eq.(lO), for a fixed value 
of J. The curve shows the logarithm of the results corresponding to the case 
Hj, for values of J < 1] = 9, and for T = 0.2 GeV and Vi = GeV. 
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Figure 9: The temperature dependence of the internal energy, as obtained 
from the calculations performed with the Hamiltonian Hj. Dashed-lines 
indicate the exact results corresponding to the unperturbed (Vi = 0) case, 
small-dashed-lines show the exact results for Vi — 0.04 GeV, and solid-hne 
shows the results of the approximations described in section III of the text, 
see Eq. (29). The validity of the approximations (solid-line curve) is limited 
to T > 0.2 GeV, as discussed in the text. The chemical potential has the 
value II = 0. 
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Figure 10: The specific lieat, as a function of temperature, for the cases 
shown in Figure 9. Results are shown following the notation given in the 
captions to Figure 9. 
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Figure 11: Temperature as a function of tlie cliemical potential. Solid-line 
indicates the exact results corresponding to the unperturbed (Vi = 0) case, 
dashed-line shows the exact results for Vi = 0.04 GeV, and small-dashed-line 
shows the results of the approximations described in section III. 
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